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The effects of the onset of deconfinement on longitudinal and transverse flow are studied. First, 
we analyze longitudinal pion spectra from fiiab — 2A GeV to ^^nn = 200 GeV within Landau's 
hydrodynamical model and the UrQMD transport approach. From the measured data on the 
widths of the pion rapidity spectra, we extract the sound velocity in the early stage of the 
reactions. It is found that the sound velocity has a local minimum (indicating a softest point in 
the equation of state, EoS) at /Sbeam = 30A GeV. This softening of the EoS is compatible with 
the assumption of the formation of a mixed phase at the onset of deconfinement. Furthermore, 
the energy excitation function of elliptic flow (v2) from isbeam — 90A MeV to y^iNN — 200 GeV 
is explored within the UrQMD framework and discussed in the context of the available data. 
The transverse flow should also be sensitive to changes in the equation of state. Therefore, the 
underestimation of elliptic flow by the UrQMD model calculation above Ex^t, = 30A GeV might 
also be explained by assuming a phase transition from a hadron gas to the quark gluon plasma 
around this energy. This would be consistent with the model calculations, indicating a transition 
from hadronic matter to "string matter" in this energy range. 
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1. Introduction 

Over the last years, a wealth of detailed data in the 20A — 160A GeV energy regime has become 
available. The systematic study of these data revealed surprising (non-monotonous) structures in 
various observables around 30A GeV beam energy. Most notable irregular structures in that energy 
regime include, 

• the sharp maximum in the K+/7r+ ratio [|l|, ^, 

• a step in the transverse momentum excitation function (as seen through {m±) — mo ) [Q, |3|], 

• an apparent change in the pion per participant ratio |^ and 

• increased ratio fluctuations (due to missing data at low energies it is unknown if this is a local 
maximum or an ongoing increase of the fluctuations) [Ql- 

It has been speculated, that these observation hint towards the onset of deconflnement already 
at 30A GeV beam energy. Indeed, increased strangeness production [|5|] and enhanced fluctuations 
have long been predicted as a sign of QGP formation [|6|, ^, |8|, ^, IC, 11, 12] within different 



frameworks and observables. The suggestion of an enhanced strangeness to entropy ratio {^^ K/n) 
as indicator for the onset of QGP formation was especially advocated in [^J]. Also the high and 
approximately constant K"^ inverse slopes of the mj- spectra above ~ 30A GeV - the 'step' - was 
also found to be consistent with the assumption of a parton <-> hadron phase transition at low SPS 
energies [|l^, |T5|]. Surprisingly, transport simulations (supplemented by recent lattice QCD (IQCD) 
calculations) have also suggested that partonic degrees of freedom might already lead to visible 



effects at ~ 30A GeV [16, 17, 18]. Finally, the comparison of the thermodynamic parameters T 



and plb extracted from the transport models in the central overlap region [ p^ with the experimental 



systematics on chemical freeze-out configurations [20, 21, ^] in the T — /ig plane do also suggest 
that a first glimpse on a deconfined state might be possible around lOA — 30A GeV. 

In the first part of this paper, we explore wether similar irregularities are also present in the 
excitation function of longitudinal observables, namely rapidity distributions. Here we employ 
Landau's hydrodynamical model and the UrQMD transport approach. In the second part, we focus 
on the excitation function of transverse flow (V2) and discuss UrQMD results in the context of the 
available data. 



2. Longitudinal flow 



It became popular to interpret relativistic heavy ion reactions with Landau's hydrodynamical 
I 28] (for recent applications of this model to relativistic nucleus-nucleus 



model [23, M 25 



interactions see [ [29[ , [30[ , p\\ |32| , |33| , |34|]). Therefore we will use this simple hydrodynamical picture 
as a baseline for the model and data comparison. The main physics assumptions of Landau's 
picture are: The collision of two Lorentz-contracted hadrons or nuclei leads to full thermalization 
in a volume of size Vnip/ ^/s. This justifies the use of thermodynamics and establishes the system 
size and energy dependence. Usually a simple equation of state /? = e/3 is assumed. Chemical 
potentials are usually assumed to vanish. The main results derived from these assumptions are: A 
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Figure 1: Left: The root mean squai'e Oy of the rapidity distributions of negatively charged pions in central 
Pb+Pb (Au+Au) reactions as a function of the beam rapidity jp. The dotted line indicates the Landau model 
prediction with = 1 /3, while the full line shows a linear fit through the data points. Data (full symbols) 
are taken from 36 3^. The statistical errors given by the experiments are smaller than the symbol 

sizes. Systematic errors are not available. Right: Speed of sound as a function of beam energy for central 



Pb+Pb (Au+Au) reactions as extracted from the data using Eq. (2.3). The statistical errors (not shown) are 
smaller than 3%. 



universal formula for the produced entropy, determined mainly by the initial Lorentz contraction 
and Gaussian rapidity distributions, at least for newly produced particles. Under the condition that 
Cs is independent of temperature, the rapidity density is given by [^, ^5|] : 

exp -— ^ (2.1) 



dy V 2^', 



with 



a. 



where is a normalisation factor and nip is the proton mass. 

Let us now analyze the available experimental data on rapidity distributions of negatively 
charged pions in terms of the Landau model. Fig. |l| (left) shows the measured root mean square 
Oy of the rapidity distribution of negatively charged pions in central Pb+Pb (Au+Au) reactions as 
a function of the beam rapidity. The dotted line indicates the Landau model predictions with the 
commonly used constant sound velocity = 1/3. The full line shows a linear fit through the data 
points, while the data points ||3, 32, ^ are depicted by full symbols. 



At a first glance the energy dependence looks structureless. The data seem to follow a linear 
dependence on the beam rapidity yp without any irregularities. However, the general trend of 
the rapidity widths is also well reproduced by Landau's model with an equation of state with a 
fixed speed of sound. Nevertheless, there seem to be systematic deviations. At low AGS energies 



and at RHIC, the experimental points are generally underpredicted by Eq. (2.2), while in the SPS 
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Figure 2: The root mean square of the rapid- 
ity distribution of negatively charged Pions in cen- 
tral AuH-Au/PbH-Pb and Proton-nProton reactions as a 
function of the center of mass energy. UrQMD calcu- 
lations for AuH-Au/PbH-Pb are denoted by full circles, 
the pp results are shown by open squares. The pre- 
diction from Landau's model is given by the line (Eq. 



2.1). Data |32| are depicted by full diamonds. 



energy regime Landau's model overpredicts the widths of the rapidity distributions. Exactly these 
deviations from the simple Landau picture do allow to gain information on the equation of state of 



the matter produced in the early stage of the reaction. By inverting Eq. ( |2.2| ) we can express the 
speed of sound in the medium as a function of the measured width of the rapidity distribution: 



4ln{./s^/2mp 



■ + 



\ 



3 a2 



+ 1 



(2.3) 



Let us now investigate the energy dependence of the sound velocities extracted from the data. 
Fig. H (right) shows the speed of sound as a function of beam energy for central Pb-i-Pb (Au-i-Au) 



reactions as obtained from the data using Eq. ( |2.3| ). The sound velocities exhibit a clear mini- 
mum (usually called the softest point) around a beam energy of 30A GeV. A localized softening of 
the equation of state is a long predicted signal for the mixed phase at the transition energy from 
hadronic to partonic matter [ p8| , p9| , pO| ]. 

To test this hypothesis the same observable has also been calculated using the UrQMD model 



(v2.2) [4^1, 42]. This transport model takes into account the formation and multiple rescattering 
of hadrons and dynamically describes the generation of pressure in the hadronic expansion phase. 
It involves also interactions of (di-)quarks, however gluonic degrees of freedom are not treated 
explicitly, but are implicitly present in strings. This simplified treatment is generally accepted to 
describe Proton-Proton and Proton-nucleus interactions. 

As depicted in Fig. ^ the UrQMD predictions (full circles) for the rapidity widths of nega- 



tively charged pions in Au-i-Au (Pb-i-Pb) reactions are in line with the experimental data ||32|] (full 
diamonds) and Landau's hydrodynamical model (full line). A rather surprising observation is that 
the calculated rapidity widths of in pp interactions (open squares) are identical to the AA re- 
sults. Together with the previous discussion, it seems that the equation of state in the transport 
model is also soft in the SPS regime. Thus, the nature of the softest point remains unclear. 
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Figure 3: Left: The calculated energy excitation function of elliptic flow of charged particles in 
Au+Au/Pb+Pb collisions in mid-central collisions (b=5-9 fm)with \y\ < 0.1 (black line). This curve is com- 
pared to data from different experiments for mid-central collisions. For E895 ^vj] and FOPI [EsJ there is 
the elliptic flow of protons and for NA49 [Hi it is the elliptic flow of pions. For E877, CERES H' 
PHENIX @, PHOBOS [|^] and STAR there is data for the charged particle flow. The dotted line in 
the low energy regime depict UrQMD calculations with the mean field |]60[|. Right: Calculated fraction of 
energy density in unformed hadrons with \y\ < 0.5 and in a cylindrical volume with transverse radius r — 3 
fm and length h^3 /Ycm fm as a function of the beam energy for central Pb-nPb (AuH-Au) reactions. 



3. Transverse flow 

Let us now look at the dynamics of the system perpendicular to the beam direction. The 
transverse flow is intimately connected to the pressure gradients. Therefore, it is sensitive to the 
equation of state (EoS) and might be used to search for abnormal matter states and phase transitions 



[ ]43[ , [44[ , ft5[ ]. Especially the second coefficient of the Fourier expansion of the azimuthal distribution 
of the emitted particles (V2) is a valuable tool to gain insight into the expanding stage of the fireball. 

The excitation function of charged particle elliptic flow is compared to data over a wide energy 
range (Fig. |3|(left)), i.e from ^beam = 90A MeV to ^ysNN = 200 GeV. The squeeze-out effect at low 
energies and the change to in-plane emission at higher energies is nicely observed in the excitation 
function. The symbols indicate the data for charged particles from different experiments. Note 
however, that in the low energy regime there are only experimental data points for protons. For 
beam energies below 2A GeV most of the charged particles are also protons because there is not 
enough energy to produce many new particles. Going to higher energies the elliptic flow of pions 
and charged particles are very similar. The rapidity cut of \y\ < 0.1 has been used for the whole 
energy range despite the fact that the data at higher energies is within \y\ < 0.5. This has been 
done to avoid too much changes in the parameters and this choice gives reasonable results over the 
whole energy range. We have checked that the results at higher energies are not affected by the 
choice of this narrower rapidity window. 

At low energies ^beam ~ 0.1 — 6A GeV the squeeze-out effect, i.e. the elliptic flow out-of- 
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plane, is clearly seen in the data as well as in the calculations, especially when the mean field is 
considered. At such energies, it is well known that both the mean field and the two-body collision 
are equally important to reproduce quantitatively the experimental results [56, 57, 58 1. In this paper 
we adopt a hard equation of state with momentum dependence (HMw) which was updated recently 
in the UrQMD model [|9|, |6^]. 

In the SPS regime the model calculations are quite in line with the data, especially with the 
NA49 results. For a more detailed discussion of directed and elliptic flow results from UrQMD-2.2 



the reader is referred to \ p\\ p2\ p3| ]. Above ^lab = 160A GeV the calculation underestimates the 
elliptic flow. At the highest RHIC energy there are about 5% flow in the data while the model 
calculation provides only half of this value. This can be explained by assuming a lack of pressure 
in the transport model at these energies. 

It is possible that above the energy range about Eiab = 30 AGeV partonic interactions have to 
be taken into account to describe the data as suggested in [18, 64, 65]. How can we analyse this 
question in the model, since there are no partonic degrees of freedom explicitly incorporated? In the 
current model exists a formation time for hadrons produced in the string fragmentation. The leading 
hadrons of the fragmenting strings contain the valence quarks of the original excited hadron. These 
(di-)quark string ends are allowed to interact during their formation time with a reduced cross 
section defined by the additive quark model. Other "pre-hadrons" from the fragmenting string are 
not allowed to interact before the coalescence of the produced quarks. Thus, because the unformed 
hadrons do not interact with others during their formation time, the effective pressure is reduced 
and only build up from the density of the formed hadrons. 

To illuminate this, we have calculated the energy density during heavy ion collisions at dif- 
ferent beam energies. From this, we extract the time corresponding to the maximum value of the 
total energy density. Fig. |3| (right) shows the fraction of the energy density that is deposited in the 
"unformed hadrons" {epartonic / {^panomc + ^hadmnic))- I-C- all string fragments within their formation 
time are dubbed as "partonic". The fraction of Epartonk starts at zero for low energies and then rises 
fast to almost 100 %. Note that this fraction reaches 90 % already around 30 AGeV beam energy, 
similar to the energy region where a phase transition is expected. As one can see, the energy den- 
sity of the formed hadrons {Ehadnmic) is much smaller than the total value, therefore the effective 
pressure of the formed hadrons alone in the model seems to be too small to generate enough V2. 
Thus, this finding supports the interpretation of the need for initial pressure from non-hadronic 
matter already at low SPS energies. 



4. Conclusion 

In conclusion, we have explored the excitation functions of the rapidity widths and of elliptic 
flow of (negatively) charged pions in Pb-i-Pb (Au-i-Au) collisions. The following observations can 
be made: 

• The rapidity spectra of pions produced in central nucleus-nucleus reactions at all investigated 
energies can be well described by single Gaussians. 

• The energy dependence of the width of the pion rapidity distribution follows the prediction 
of Landau's hydrodynamical model if a variation of the sound velocity is taken into account. 
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• The speed of sound excitation function exti-acted from the data has a pronounced minimum 
(softest point) at ^beam = 30A GeV. 

• The UrQMD model describes the rapidity widths data well, but underestimates the elliptic 
flow from the higher SPS energy on. 

• The softest point coincides with the rapid rise of "partonic" degrees of freedom in the present 
model. 

• This softest point might be due to the formation of a mixed phase indicating the onset of 
deconfinement at this energy. 

Further explorations of this energy domain is needed and can be done at the future FAIR 
facihty and by CERN-SPS and BNL-RHIC experiments. 

Acknowledgements 

This work was supported by GSl and BMBF This work used computational resources pro- 
vided by the Center for Scientific Computing at Frankfurt (CSC). 

References 

[1] S. V. Afanasiev et al. [The NA49 Collaboration], Phys. Rev. C 66 (2002) 054902 
[arXiv:nucl-ex/0205002] . 

L2J M. Gazdzicki et al [NA49 Collaboration], J. Phys. G 30 (2004) S701 [arXiv:nucl-ex/0403023]. 

[3] C. Blume, J. Phys. G: Nucl. Part. Phys. 31, S57 (2005) 

[4] C. Roland [NA49 Collaboration], J. Phys. G 31 (2005) S1075. 

[5] P Koch, B. Muller and J. Rafelski, Phys. Rept. 142 (1986) 167. 

[6] M. Bleicher, S. Jeon and V. Koch, Phys. Rev C 62 (2000) 061902 [arXiv:hep-ph/0006201]. 

[7] S. Jeon, L. Shi and M. Bleicher, Phys. Rev. C 73, 014905 (2006) [arXiv:nucl-th/0506025]. 

[8] E. V. Shuryak and M. A. Stephanov, Phys. Rev C 63 (2001) 064903 [arXiv:hep-ph/0010100]. 

[9] H. Heiselberg and A. D. Jackson, Phys. Rev C 63 (2001) 064904 [arXiv:nucl-th/0006021]. 

[10] B. Muller, Nucl. Phys. A 702 (2002) 281 [arXiv:nucl-th/01 11008]. 

[11] M. Gazdzicki, M. I. Gorenstein and S. Mrowczynski, Phys. Lett. B 585, 1 15 (2004) 

[arXiv:hep-ph/0304052] . 

[12] M. I. Gorenstein, M. Gazdzicki and O. S. Zozulya, Phys. Lett. B 585, 237 (2004) 
[arXiv:hep-ph/0309142]. 

[13] M. Gazdzicki and M. I. Gorenstein, Acta Phys. Polon. B 30, 2705 (1999). 

[14] M. 1. Gorenstein, M. Gazdzicki and K. A. Bugaev Phys. Lett. B 567, 175 (2003) 
[arXiv:hep-ph/0303041]. 

[15] Y. Hama, F. Grassi, O. Socolowski, T. Kodama, M. Gazdzicki and M. Gorenstein, Acta Phys. Polon. 
B 35 (2004) 179. 



7 



Longitudinal flow and onset of de confinement 



Marcus Bleicher 



[16] H. Weber, C. Ernst, M. Bleicher et al, Phys. Lett. B 442, 443 (1998). 
[17] E. L. Bratkovskaya etal., Phys. Rev. Lett. 92, 032302 (2004) 
[18] E. L. Bratkovskaya etal., Phys. Rev. C 69, 054907 (2004) 

[19J L. V. Bravina et al, Phys. Rev. C 60, 024904 (1999), Nucl. Phys. A 698, 383 (2002). 

[20] P Braun-Munzinger and J. Stachel, Nucl. Phys. A 606, 320 (1996). 

[21] R Braun-Munzinger and J. Stachel, Nucl. Phys. A 638, 3 (1998). 

[22] J. Cleymans and K. Redhch, Phys. Rev. C 60, 054908 (1999). 

|23J E. Fermi, Prog. Theor. Phys. 5, 570 (1950). 

[24] L. D. Landau, Izv. Akad. Nauk Ser Fiz. 17, 51 (1953). 

[25] S. Z. Belenkij and L. D. Landau, Usp. Fiz. Nauk 56, 309 (1955). 

[26] P Carruthers and M. Doung-van, Phys. Rev. D 8, 859 (1973). 

[27] P Carruthers, Annals N.Y.Acad.Sci. 229, 91 (1974). 

[28] P A. Carruthers, LA-UR-81-2221 

[29] J. Stachel and P Braun-Munzinger, Phys. Lett. B 216, 1 (1989). 
[30] P Steinberg, arXiv:nucl-ex/0405022. 
[31] M. Murray, arXiv:nucl-ex/0404007. 

[32] G. Roland, Talk presented at Quark Matter 2004, see proceedings. 
[33] M. Bleicher, arXiv:hep-ph/0504207. 
[34] M. Bleicher, arXiv:hep-ph/0509314. 
[35] E. V. Shuryak, Yad. Fiz. 16, 395 (1972). 

[36] J.Klay et al. [E895 Collaboration], Phys. Rev. C 68, 054905 (2003) 

[37] l.G. Bearden et al. [Brahms Collaboration], Phys. Rev. Lett. 94, 162301 (2005) 

[38] C. M. Hung and E. V. Shuryak, Phys. Rev. Lett. 75, 4003 (1995) [arXiv:hep-ph/9412360]. 

[39] D. H. Rischke, Y. Pursun, J. A. Maruhn, H. Stoecker and W. Greiner, Heavy Ion Phys. 1, 309 (1995) 
[arXiv:nucl-th/9505014]. 

[40] J. Brachmaim, A. Dumitru, H. Stoecker and W. Greiner, Eur. Phys. J. A 8, 549 (2000) 
[arXiv:nucl-th/9912014]. 

[41] S.A. Bass et al, Prog. Part. Nucl. Phys. 42, 255 (1998). 

[42] M. Bleicher et al., J. Phys. G 25, 1859 (1999). 

[43] H. Stocker, J. A. Maruhn and W. Greiner, Z. Phys. A 290 (1979) 297. 

[44] J. Hofmann, H. Stocker, U. W. Heinz, W. Scheid and W. Greiner, Phys. Rev. Lett. 36 (1976) 88. 
[45] H. Stocker and W. Greiner, Phys. Rept. 137 (1986) 277. 

[46] C. Pinkenburg et al. [E895 Collaboration], Prepared for Centennial Celebration and Meeting of the 
American Physical Society ( Combining Annual APS General Meeting and the Joint Meeting of the 
APS and the AAPT), Atlanta, Georgia, 20-26 Mar 1999 



8 



Longitudinal flow and onset of de confinement 



Marcus Bleicher 



[47] P. Chung etal. [E895 CoUaboration], Phys. Rev. C 66 (2002) 021901 [arXiv:nucl-ex/01 12002]. 
[48] A. Andronic etal. [FOPI CoUaboration], Phys. Lett. B 612 (2005) 173 [arXiv:nucl-ex/0411024]. 
[49] C. Alt etal. [NA49 Collaboration], Phys. Rev. C 68 (2003) 034903 [arXiv:nucl-ex/0303001]. 
L50J K. FiUmonov et al. [CERES/NA45 Collaboration], arXiv:nucl-ex/0109017. 

[51] J. Slivova [CERES/NA45 Collaboration], Nucl. Phys. A 715 (2003) 615 [arXiv:nucl-ex/02120131. 

[52] S.l.Esumi, J. Slivova, J. Milosevic for CERES Collaboration SFIN, year XV, Series A: Conferences, 
No. A2(2002) 

[53] S. Esumi [PHENIX Collaboration], Nucl. Phys. A 715 (2003) 599 [arXiv:nucl-ex/0210012]. 

[54] S. Manly et al. [PHOBOS Collaboration], Nucl. Phys. A 715 (2003) 611 [arXiv:nucl-ex/0210036]. 

[55] R. L. Ray [STAR Collaboration], Nucl. Phys. A 715 (2003) 45 [arXiv:nucl-ex/021 1030]. 

[56] R Danielewicz, Nucl. Phys. A 661 (1999) 82 [arXiv:nucl-th/9907098]. 

[57] P. Danielewicz, R. A. Lacey, P. B. Gossiaux, C. Pinkenburg, P. Chung, J. M. Alexander and 
R. L. McGrath, Phys. Rev. Lett. 81 (1998) 2438 [arXiv:nucl-th/9803047]. 

[58] Q. b. Pan and P Danielewicz, Phys. Rev. Lett. 70 (1993) 2062 [Erratum-ibid. 70 (1993) 3523]. 

[59] Q. Li, Z. Li, S. Soff, M. Bleicher and H. Stocker, J. Phys. G: Nucl. Part. Phys. 32 (2006) 151 

[60] Q. f. Li, Z. X. Li, S. Soff, M. Bleicher and H. Stocker, J. Phys. G 32 (2006) 407 
[arXiv:nucl-th/0601047]. 

[61] H. Petersen, Q. Li, X. Zhu and M. Bleicher, Phys. Rev. C in print [arXiv:hep-ph/0608189]. 

[62] X. 1. Zhu, M. Bleicher and H. Stoecker, J. Phys. G 32 (2006) 2181 [arXiv:nucl-th/0601049]. 

[63] X. 1. Zhu, M. Bleicher and H. Stoecker, Phys. Rev. C 72, 06491 1 (2005) [arXiv:nucl-th/0509081]. 

[64] H. Weber et al, Phys. Lett. B 442 (1998) 443 [arXiv:nucl-th/9808021]. 

[65] X. Zhu, H. Petersen and M. Bleicher, AlP Conf. Proc. 828 (2006) 17. 



9 



